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1 Introduction 



The main problem of stationary quantum chaos is the following: Given a quantum 

Hamiltonian (operator) H with infinitely many hound states, being the quantized 
object based on its classical analog H = H{q,p), the Hamiltonian H as a function 
of the N coordinates q and momenta p, what are the geometrical and statistical 
properties of the eigenf unctions, their Wigner functions, of the energy spectra and 
of the matrix elements of other quantal observables. 

Of course, we have the (stationary) Schrodinger equation, which sometimes we 
can even solve analytically, e.g. in cases of analytically solvable one- dimensional 
potentials or in cases of separable A^-dimensional potentials (Landau and Lifshitz 
1997), and sometimes we can solve the underlying eigenvalue problem numerically, 
at least in principle. However, the analytically solvable problems are very untypical, 
although quite important, because we can use them to explore their neihbourhood 
(in the functional space of Hamiltonians H) by means of a large variety of per- 
turbational techniques, and such a neighbourhood includes classically nonintegrable 
systems which are typical (generic) . This is in analogy of the KAM systems in classi- 
cal mechanics (Gutzwiller 1990). The numerical techniques can be applied to almost 
all systems, but it turns out that as soon as the system is not (classically) integrable 
and solvable, also the numerical techniques must be quite sophisticated, especially 
if we ask for high lying eigenstates. To know the Schrodinger equation and to have 
the potentiality of solving it helps as httle as the analogous situation in the classical 
dynamics where the potentiality of solving the Hamilton- Jacobi equation does not 
help very much in studying the global, qualitative and quantitative, properties of 
motion in generic nonintegrable Hamiltonian systems. This has been realized for the 
first time by Henri Poincare, who has shown (see Poincare 1993, Goroff 1993) that 
the gravitational three-body problem is indeed nonintegrable and this broken inte- 
grability can no longer warrant the existence of invariant tori everywhere in phase 
space, thereby giving the way to true chaotic motion, which cannot be embedded 
into smooth A^-dimensional invariant surfaces. The emergence of classically chaotic 
motion gives rise to the notion of qualitative dynamics, which is apt to embrace the 
richness and variety of chaotic behaviour. To achieve that new methods are needed, 
both analytically and numerically. Concepts of Surface-of-Section (SOS) and similar 
ones become indispensible to work in the new science of nonlinear dynamics, dealing 
with the dynamical systems described by the set of M ordinary first order differ- 
ential equations in their phase space. For Hamiltonian systems we have M — 2N 
Hamilton equations. On the quantum side we face a precisely analogous problem: 
Given the (stationary) Schrodinger equation of a A^-dimensional quantum system, 
whose classical analogue is not only nonseparable but also nonintegrable and thus 
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chaotic, we can hardly see the structure of the solutions (eigenstates, described by 
the wave function and the corresponding Wigner function) and their global proper- 
ties. New approach is necessary, including the numerical one, to see and classify all 
the possible types of behaviour. On the theoretical analytical side the semiclassical 
methods are quite essential, and this line of thoughts goes back to the pioneering 
and classical work of Gutzwiller (Gutzwiller 1990 and references therein) and Per- 
cival (1973), further developed by many workers in classical and quantum chaos 
(Chirikov 1979, Casati and Chirikov 1994, Berry 1983, Giannoni et al 1991, Haake 
1991, Bohigas and Giannoni 1984, Bohigas 1991). For a recent excellent review, 
covering not only quantum chaos, but also all related theoretical and experimental 
branches of physics see the paper by Weidenmiiller and coworkers (Guhr et al 1998). 

The purpose of this paper is to review the main methods and results in the field of 
quantum chaos, i.e. the study of the solutions of the Schrodinger equation connected 
with the classically nonintegrable and chaotic systems. 



2 The main assertion of stationary quantum chaos 

The main assertion of stationary quantum chaos is the following answer to the main 
problem of quantum chaos in the semiclassical limit of sufficiently small h: 

2.1 Classical integrability 

The case (I) of classically integrable quantal systems H: 

If H is classically integrable, then the wave function is locally a superposition of 
a finite number of plane waves, the number of directions of the wave vector being 
equal to the number of possible momentum vectors p through the coordinate point 
q. Globally, the probability density is equal to the classical probability density 
obtained by projecting the A^-dimensional invariant torus onto the configuration 
space q, up to within the resolution scale of the order of one de Broglie wave length. 
The corresponding Wigner function^ (q, p) of the eigenstate is a delta function on 
the invariant torus labeled by the quantized classical action variable I(q, p) = In, 
where n is the quantum number multi-index n = (rii, n2, . . . , n^r) denoting the 
Maslov (EBK) quantized invariant A^-torus, so that 

^For the definition and properties of Wigner functions see section 3. 
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W{ci,p) = ^SM^,p) -1^) (1) 

where 6f is the /-dimensional Dirac deha function, and in our case f = N. The 
eigenvalues of H, i.e. the eigenenergies, in a small interval, after unfolding, that is 
after reducing the mean energy level spacing to unity, obey (typically) the Poissonian 
statistics (Robnik and Veble 1998): The probability E{k,L) of observing k levels 
inside an interval of lentgh L is given by^ 



Ik 

^integrable ~ ^Poissonianiki = -^j-exp( Z/) (2) 

The untypical cases have measure zero, and are characterized by some number the- 
oretic special properties like e.g. the rectangle billiards with rational squared sides 
ratio. 



2.2 Classical ergodicity 



The case (E) of classically ergodic quantal systems H: 

If H is classically ergodic system, then the wave function is locally a superposition 
of infinitely many plane waves, the directions of the wave vector k being isotrop- 
ically distributed on a A^-dimensional sphere of radius ^2m{E — y(q)), if the 
Hamiltonian is iJ(q, p) = p^/(2m) + ^(q), where m is the mass and V the poten- 
tial energy. Due to the ergodicity the phases of the plane waves are assumed to be 
random (Berry 1977a,b), which implies that the wave amplitude V'(q) is a Gaussian 
random function. The random phase assumption, however, can break down in vicin- 
ity of isolated unstable classical periodic orbits or families of such orbits, where we 
observe the scars (Heller 1984), i.e. the regions of enhanced probability density 
Thus, to the leading approximation we have the microcanonical Wigner function of 
almost a// eigenstates (Shnirelman 1979, Berry 1977a,b, Voros 1979), 



I^1q.p)- ,,/\lf-^''^'^'' „ (3) 

j dNqd^pdi{E - iZ(q,p)) 



■^It can be shown that knowledge of all E(fc, L)-statistics is equivalent to the complete knowledge 
of all n-point correlation functions. The calculation of E(A;,L), however, is much easier than the 
calculation of correlation functions etc., since there is no binning and other advantages. This has 
been pointed out very clearly by Aurich et al (1997). See also the book by Mehta (1991). 
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The eigenenergies of -ff in a small interval, after unfolding, obey the predictions of 
classical Random Matrix Theories (RMT), namely the statistics of the eigenvalues 
of the ensemble of orthogonal Gaussian matrices (GOE) or of unitary Gaussian 
ensembles (GUE) (depending on the existence or nonexistence of an antiunitary 
symmetry), introduced by Wigner, and also of the COE/CUE of Dyson (see e.g. 
the book by Haake (1991)). (We consider only systems with classical analog H 
and therefore ignore spin and GSE). This assertion has been proposed by Bohigas, 
Giannoni and Schmit (1984). It implies that E(k, L) statistics must obey the RMT 
laws, the so-called BGS-Conjecture 

Eergodic{k, L) = ERMT{k, L) (4) 
2.3 Classically mixed systems 

The case (M) of classically mixed (generic) quantal systems H: 

If H is classically mixed system, then we can distinguish between regular and irregu- 
lar states. Percival (1973) was the first to propose such a qualitative characterization 
of eigenstates. The regular states are associated with classical invariant tori (semi- 
classically EBK/Maslov quantized tori, to the leading semiclassical approximation), 
and the chaotic states are associated with chaotic components. This view has been 
made more quantitative in the work of Berry and Robnik (1984). The Berry- Robnik 
picture rests upon the The Principle of Uniform Semiclassical Condensation (PUSC, 
see section 3), which states that the Wigner functions of quantal states in the limit 
h ^ Q become positive definite, and since they are mutually orthogonal, they must 
"live" on disjoint supports, and the phase space volume (Liouville measure) of each 
of them is of the order of {2tt%)^ . See section 3 and e.g. (Robnik 1997). The 
question is, what is the geometry of the object on which they "condense", and the 
answer - as a conjecture - is: uniformly on a classical invariant object (Berry 1977a,b, 
Robnik 1988, 1995, 1997). Therefore we have regular and irregular states. The as- 
sumption is that there is no correlation between the spectral sequences (regular and 
a series of irregular states). If > 3 we have only one chaotic component (the 
Arnold web of chaotic motion pervades the entire phase space - energy surface - and 
is dense, i.e. its closure is the energy surface) and one associated irregular sequence 
of eigenstates, whereas in = 2 we have many, even infinite number of sequences 
of irregular states, of smaller and smaller invariant measure, each sequence being 
associated with one chaotic component. It is thus assumed that the Wigner function 
of a regular state is of type (|l|), whilst for irregular states, and generally, it is 
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Wia v) = ^/(F(q,p))x..(q,p) , . 

where Xuji'i, p) is the characteristic function of the invariant component, labeled by 
u, being a (either smooth or nonsmooth, generally possibly also fractal) subset of 
the smooth (2A^ — /)-dimensional invariant surface defined by the / implicit equa- 
tions (global integrals of motion), namely F(q, p) = 0, where F = {Fi, F2, . . . , Fj). 
The characteristic function Xu]{^.iV) is defined to have value unity on uj and zero 
elsewhere. The integer number / can be anything between 1 (ergodic system) and 
(integrable system). 

Obviously, the formula (^ is the most general expression for a condensed Wigner 
function of a (pure) eigenstate. It generalizes the cases (I) and (E). Namely, if we 
have ergodicity, then / = 1, we put Fi(q, p) = £^ — -ff(q, p), and u = entire energy 
surface, and we recover equation (^). In the other extreme (I), we have N global 
integrals of motion in involution, and so F(q, p) = — I(q, p), and u = the invari- 
ant torus labeled by !„, and we recover equation (|l]) of case (I). In the most general 
case, therefore, formula (|) applies. Obviously, W is normalized (see next section), 

J d^cid^pW{ci, p) = 1 (6) 

For generic (mixed) systems the most typical case is / = 1, -Fi(q, p) = E — iJ(q, p) 
and UJ is a (nonsmooth, typically fractal, chaotic) subset of the energy surface Fi. 
We write down this most important case explicitly: 

W(q p) = S^jE - H{ci,p))xM,p) 

^' J rf^qd^p5i(E - Hiq, p))x.(q, p) ^ ^ 

It is important to know the relative invariant (Liouville) measure of chaotic and 
regular eigenstates because the Hilbert space of a mixed Hamiltonian system is split 
into regular and irregular eigenstates, in the strict semiclassical limit, precisely in 
proportion to the classical invariant measure of the integrable component (invariant 
tori) and of the irregular components. 

It is quite obvious by looking at the equation (|^) that the invariant Liouville measure 
of a subset u of the energy surface is equal to 

. _ Jd^qd^p5,{E-H{q,p))xM.P) .o^ 
^^""^ Jd^qd^p6^{E-H{q,p)) 



6 



The relative invariant Liouville measure of the regular components will be denoted 
by pi, and the measures of chaotic components (ordered in sequence of decreasing 
measure) by p2, Ps, ■ ■ ■ , Pm, where m = cxd for = 2 and m = 2 for iV > 3, as 
already explained. In section 4 I shall explain how one can calculate the measures 
P2,P3, • ■ ■• 

Assuming the above mentioned absence of correlations pairwise between m spec- 
tral sequences, due to the fact that they have disjoint supports and thus do not 
interact, where m is infinite for N = 2 and 2 for > 3, the spectral statistics can 
be written as 

m 

ki+k2 + ...+km=k j=l 

which is a manifestation of Berry- Robnik (1984) picture. Here Ej(fc, L) is Epoisson{k, L) 
for j = 1, and EptMrik, L) for j = 2, 3, . . . , m. See cases (I) and (E), equations (Q) 
and d^). The picture is based on the reasonable assumption that (after unfolding) 
the mean density of levels in the j-th sequence of levels is pj, simply applying the 
Thomas- Fermi rule of filling the phase space volume with elementary cells of size 
in the thin energy shell embedding the corresponding subset oj. Therefore, 
please note that the second argument of Ej(/c, L) is weighted precisely by the classi- 
cal relative invariant measure of the underlying invariant component. Also, if there 
were several regular (Poissonian) sequences they can be lumped together into a sin- 
gle Poissonian sequence (which we traditionally label by 1 with relative invariant 
measure pi): It is easy to show, that if ai, 02, are positive real numbers and 
[3 being their sum, (3 = ai + a2 + ■ ■ ■ + on, then for all fc, and L, 

Epoisson{k, PL) = ^ Epoisson{ki, aiL)Epoisson{k2, a2L) . . . Epoisson{ki, aiL) 

ki+k2 + ...+ki=k 

(10) 

by simply using the definition of Epdssonik, L) of equation (|^). Thus, we have some 
kind of a central limit theorem, saying that the statistically independent superposi- 
tion of Poisson sequences results in a Poisson sequence, such that the total density of 
Poissonian levels (5 is equal to the sum of the partial level densities Oj, j = 1, 2, . . . , / . 

The case (M) is the most general one, and as the limiting extreme cases includes 
cases (I) and (E). 
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2.4 Limitations of the universality 



There are two important limitations of the above stated asymptotic behaviour as 
^ — >■ 0, when h is not yet small enough: One is the existence of the outer energy 
scale, and the other one is the localization phenomena. As for the first, it has been 
shown by Berry (1985), applying the semiclassical Gutzwiller periodic orbit theory 
(Gutzwiller 1990 and the references therein), that at energy scales (after unfolding!) 
L > Lmax we do not have the universality but typically a saturation, i.e. E(fc, L) 
statistics at L larger than 



l^max ^^^^^^ liij 

where ( Ai?) is the mean energy level spacing, and Tq the period of the shortest clas- 
sical periodic orbit in the dynamical system if(q, p), deviate from their universal 
behaviour of cases (I) and (E), equations (Q) and (^. Instead, e.g. the sigma and 
delta statistics become constant. (For definition and inter-relationship see section 
3.) However, please observe that as h goes to zero, also Lmax goes to infinity as 
a power h^'^^^, so that in the semiclassical limit the universality region becomes 
larger oc fi~^^^. 

As for the second limitation bordering the universal behaviour we comment the 
following. If the value of the (effective) Planck constant h is not sufficiently small, 
then the eigenf unctions might not be fully extended in the sense of the correspond- 
ing Wigner functions obeying the equations (|I]), (^ and (^, but can be localized 
(not uniformly extended) on the classical invariant object on which they condense. 
Such a deviation from the ultimate limiting semiclassical behaviour is therefore a 
manifestation of the localization phenomena in stationary eigenstates of autonomous 
Hamiltonian systems, and is manifested also in the spectral statistics. For example, 
if we have a classically ergodic system, but with very slow chaos (very large diffu- 
sion time), we shall observe strongly localized states, mimicking a regular integrable 
system. In such an extreme case of localization we shall observe Poissonian spec- 
tral statistics rather than GOE/GUE. Depending on the strength of localization we 
shall therefore be able to see transition from Poissonian to GOE/GUE behaviour 
in an ergodic ssystem. In the intermediate crossover regime we observe Brody-like 
behavior with fractional power law level repulsion. In a KAM regime the same is 
true if the effective h is not small enough to resolve the structure of small chaotic 
components. We shall describe these phenomena in the subsequent sections. 

However, just briefly, a qualitative comment is in order at this place. The rele- 
vant criterion for localization is, that the so called break time or Heisenberg time, 
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defined through 



tbreak - tn - (^2) 

where (A£^) is the mean level spacing, must be shorter than the diffusion time, 
so then we have strongly localized states, whilst in the opposite extreme we have 
strongly extended states. The reason is very simple: Quite generally quantum me- 
chanics (of a suitably chosen initial wave packet) follows classical dynamics (of a 
suitably chosen ensemble of initial conditions) up to the break time, after which 
the interference phenomena set in, resulting typically in destructive interference, 
and thus in the stop of diffusion, which means localization (before the entire phase 
space has been conquered). For example, in two-dimensional billiards, (AE) is just 
constant, so the break time is constant and independent of energy, whilst the clas- 
sical diffusion time, even if very large at small energies E, decreases with energy 
as const./ y/E, so ultimately, as — > cxo, we shall find the extended states and 
then the general picture of case (M) is applicable, which, of course, as the extreme 
cases, includes (I) and (E). However, the phenomena of localization, including the 
scars, are extremely important and as we have sketched above, they are related to 
the important time scales which control the finite time structure and behaviour of 
classical dynamics, especially the transport times and so on. 



2.5 Distribution and fluctuation properties of transition prob- 
abilities 

Finally, in this subsection we should comment on the statistical properties of the 
matrix elements of other observables in the eigenbasis of an integrable, ergodic and 

mixed system. The main work in this direction has been done by Feingold and Peres 
(1986) for the ergodic case, and this has been been generalized to integrable and 
mixed systems by Prosen and Robnik (1993a). 

The expectation values and generally the matrix elements of other reasonable observ- 
ables (Hermitian operators having a classical limit) have been little studied (Fein- 
gold and Peres 1986, Alhassid and Levine 1986, Wilkinson 1987, 1988). One well 
known result concerns the fluctuation properties of generalized intensities (squares 
of matrix elements) within the framework of random matrix theories, namely the 
Porter-Thomas distribution (Brody et al 1981), which has been experimentally ob- 
served and suggested by Porter and Thomas (1956) in the context of nuclear physics. 
We expect that this fluctuation law applies also in classically ergodic systems with 
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few freedoms. The main motivation of our work (Prosen and Robnik 1993a) was to 
explain this and to find the appropriate generahzation for Hamiltonian systems in 
the transition region of mixed dynamics. 

In order to study the fluctuation properties of generahzed intensities one must be 
able to clearly separate the smooth mean part of the intensities as the function of 
frequency (= energy difference between the final and initial state/h) from its fluc- 
tuating part. So, given the frequency of the intensity we ask what is its mean value 
and which is the distribution of its fluctuating part in units of the mean value. In the 
classically ergodic case Feingold and Peres (1986) propose a formula expressing the 
mean intensities in terms of the power spectrum of the given observable taken over 
a dense chaotic classical orbit. In deriving this result they rely on the Shnirelman 
theorem (1979) expressing the quantum expectation value of a reasonable operator 
as the classical microcanonical average. This theorem is obvious once one has in 
mind that the Wigner distributions of the eigenstates of a classically ergodic system 
in the semiclassical limit are just microcanonical distributions (Berry 1977a,b, Voros 
1979), equation (|^). In order to rederive Feingold- Peres formula and to generalize 
it we first point out that the Shnirelman theorem applies also to the states in the 
regular and mixed regime if the classical average is taken over the relevant classical 
invariant ergodic component, which supports the corresponding semiclassical eigen- 
state. This can be an invariant torus, a chaotic component, or the entire energy 
surface. 

Following Feingold and Peres (1986) we start by looking at the following sum over 
eigenstates k of eigenenergies for the transition elements Aj^ = {j\A\k) 



Y^exp{z{E,-E,)t/fi)\A,,\' = 5:(j|e^^^*/^i|A;)(fc|e-^^^-*/'^i|j)= (13) 

k k 

= (j|e^^*/^ie-^*/^i|j) = (j|i(t)i(0)|j) 



Now we apply the generalized Shnirelman theorem, stating that in the semiclassical 
limit this is equal to the classical average 



over the invariant ergodic component labeled by j which supports the semiclassical 
state Using the ergodicity on the given invariant component this two-point 
autocorrelation function can be expressed as the time average along a classical dense 
orbit (dense in the given invariant component, which e.g. can be an invariant torus, 
or a chaotic component, or the entire energy surface) 



Q(t) = {Arnm. 



3 



(14) 




r/2 



(15) 



■r/2 
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Next we replace the sum by the integral / dEkp{Ek), where p{E) is the density 
of states, and perform the Fourier transform and obtain 



,2, - Ej)/K) 



where the state j is fixed and the average is taken over states k within a thin 
energy shell of thickness of few mean level spacings. Here 



Sjico)^ I dtC{t)e'^' = \\m^ J dtA{t)e 



T/2 

-iU!t 



-T/2 



(17) 



is the power spectrum of a dense orbit in the invariant ergodic component j. HA 
has a nonvanishing mean value {A}j the S.j{uj) will have a delta spike at u; = 0, and 
this can be removed by replacing A in the above formulas hj A — {A}j. To calculate 
the actual mean values of the intensities l^jfcp we perform in the above formula (on 
the LHS) also the averaging over the j states microcanonically over the thin energy 
shell around Ej of sufficiently wide width such that the corresponding semiclassical 
states uniformly cover the energy surface, whilst on the RHS we correspondingly 
take the microcanonical average over all initial conditions j on the energy surface 
Ej. So the final formula for the mean generalized intensities is 

,^_{S{{E,-E,)/h)}E 



27rhp{E) 



By {.}e we denote the microcanonical average over the energy surface E. The ap- 
parent asymmetry in jk of this formula disappears in the semiclassical limit h ^ 0. 
In the numerical evaluations described in (Prosen and Robnik 1993a) we applied 
the above formula with {S{u)}e and p{E) being calculated on the energy surface 
placed half way between Ej and Ek, i.e. E — (Ej + Ek)/2. This choice is met to 
minimize the error at finite h. 

Knowing the average value of intensities as a function of u we can now separate 
the smooth part from its fluctuating part by renormalizing the matrix elements as 
follows 

Xjk = , • (19) 

The renormalized matrix elements Xjk are now regarded as random variable whose 
probability distribution is denoted by D{X), which by definition has unit disper- 
sion, and naturally is expected to be even function of X, D{X) = D{—X), and 
so it has zero mean. In the classically ergodic case we expect that quite gener- 
ally the matrix elements of a given operator are very well modelled by the GOE of 
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random matrix theories (Brody et al 1981) which predict Gaussian distribution for 
Dpt{X) = exp(— X^/2)/v^27r, which is equivalent to the so-called Porter-Thomas 
distribution for the intensities / = X^, namely P{I) = exp(— //2)/v27r/, see (Porter 
and Thomas 1956). In integrable cases one expects a vast abundance of at least ap- 
proximate selection rules which render most X to become zero implying that D{X) 
approaches a delta function 6{X) in the semiclassical limit. This can be seen by 
considering matrix representation of an operator in the basis of the torus quantized 
eigenstates of an integrable system, as explained in detail in (Prosen and Robnik 
1993b). In the mixed type dynamics (KAM) in the transition region between inte- 
grability and chaos we expect a continuous transition from S{X) towards Dpt{X). 
More precisely, a semiclassical formula for D{X) in such transition region has been 
derived by Prosen (1994a) by taking into account the fact that the only broadening 
of D{X) stems from the transitions between chaotic initial and chaotic final states 
belonging to the same family of the invariant ergodic components (continuously 
parametrized by the energy), while all other transitions are almost forbidden. This 
work rests upon a more detailed analysis of higher autocorrelation functions and is 
reported on in (Prosen 1994b). More details and the numerical illustration of our 
results can be found in (Prosen and Robnik 1993a) and in (Prosen 1994a,b). 



3 The Principle of Uniform Semiclassical Con- 
densation and more about the wave functions 
and statistics 

In this section we want to explain the main ingredients and arguments leading to 
the equations (|T]), (^ and (|]). The ideas go back to Berry (1977a,b), Shnirelman 
(1979), Voros (1979), Robnik (1988, 1995, 1997). To see that the quantum analogy 
of the stationary (aspects of) chaos works well it is necessary to look at the objects 
uniquely determined by given eigenstates in such a manner that one can compare 
the eigenstates to the classical states (phase portraits at given energy). This can be 
achieved by introducing the Wigner functions (transforms) of given eigenstates, e.g. 
of the wave functions. 

With this procedure we are building up a kind of the quantal phase space, in the 
spirit of the Wigner- Weyl formalism (de Groot and Suttorp 1972), in the follow- 
ing way: Let '0n(q) be the n-th wave function (eigenfunction as a solution of the 
Schrodinger problem) in the A^-dimensional configuration space with q being a po- 
sition vector, then its corresponding Wigner function (or transform) is defined as: 
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WM,P) = J c?^Xexp(-lp.X)V^„(q- |)C(q+ |) (20) 

where * denotes the complex conjugation. The Wigner function is obviously real. 
However, unlike the classical distribution functions, the Wigner functions are not 
positive definite, which is a fundamental consequence of the very nature of quantum 
mechanics: If they were positive then the quantum mechanics would be identical to 
classical mechanics. This is the essence of the Wigner theorem (about the quantal 
phase space distributions). On the other hand, the Wigner functions do have the 
correct property that they become the configurational probability density |V'(q)P 
when projected down onto the configuration space (i.e. integrating ([20| ) over the 
momenta p) and complementary, they become the momentum probability density 
|y9(p)P when integrated over the entire configuration space, as can be immediately 
verified from the definition (|20|) . Therefore Wn integrates to 1 if ipn is normalized 
to unity, i.e. 

J W„(q,p)ci^qci^p = l. (21) 
They also obey the orthogonality relation 

{27rtif J d^'qd^pWniq, p)WUq, p) = ^n™, (22) 

{6nm here is the Kronecker delta, i.e. discrete delta function, equal to 1 if n = m 
and zero otherwise) which can be understood at once by recalling that {27Th)^Wn is 
in fact the Weyl symbol of the projection operator P„ = \n >< n\. Further, it can 
be easily seen that the absolute value of Wn is bounded from above, namely (Baker 
1958, see also de Groot and Suttorp 1972, and Berry 1977b) 

|I^n(q,p)|<(^)'^ (23) 

showing that it can diverge only in the semiclassical limit when h ^ 0. One can 
also see from (p^) that (when n = m) one has 

I W^{q, p)t/^qrf^p = l/i27ihf (24) 

Therefore unlike (|2T|) the latter integral (0) can diverge as h ^ 0, and this diver- 
gence can be due to large contributions for large values of |q| and |p| or due to the 
singularities of (Berry 1977a). The latter possibility is the one that actually 
occurs, as shown by Baker (1958) 

Wn{q,p) ^ {2nhfW^{ci,p), as k ^ 0. (25) 
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Therefore in the semiclassical hmit the Wigner functions Wn become positive defi- 
nite and divergent ~ (27r^)~^, which is weak enough still to obey (|22|) with n = m. 
This orthogonality relation shows then that all pairs Wn-, Wm, with n ^ m must 
have disjoint supports. Therefore they are effectively nonzero and divergent only 
on a small piece of volume of size {2t[K)^ . In fact this semiclassical condensation 
must take place close to and on the energy shell around the classical energy surface 
of energy of total volume {27ih)^ . This is of course equivalent to the simple 
Thomas-Fermi rule of how to determine the average density of states, semiclasi- 
cally: Divide the available classical phase space volume inside the given energy En 
by {27ih)^ , and on the average this must be equal to n, which is the sequential 
quantum number, i.e. the (cumulative) number of eigenstates below the energy En- 
From these general considerations we cannot conclude more than stated. It is not 
clear a priori on what geometrical object does the Wigner function Wn condense as 
h ^ 0. Berry (1977a,b) and Voros (1979), in agreement with Shnirelman's theorem 
(1979), have suggested that in case of classical ergodicity Wn condenses uniformly 
on the energy surface, becoming the microcanonical distribution, (^. On the other 
hand, in the opposite extreme of an integrable motion, or a KAM system with in- 
variant tori, which are EBK quantized and support the quantum state \n >, Berry 
(1977a) derives from the semiclassical wave functions in coordinate space that the 
corresponding Wigner function is equal to (|I]). So here the Wigner function con- 
denses uniformly on the EBK quantized torus. 

The latter result can be easily obtained by noting that in the semiclassical limit 
the classical canonical tranformations and the quantization do commute, and thus 
can be performed directly in the space of action-angle variables, immediately yield- 
ing the above result, as shown in (Robnik 1995, Hasegawa et al 1989). 

In both extreme cases, the ergodicity (|^), and the quasi-integrability (existence of a 
quantized KAM torus) (|I|), we see that the Wigner function condenses uniformly on 
the underlying classical object, which is the invariant indecomposable component 
in the classical phase space. It seems thus very natural to elevate these findings to 
the Principle of Uniform Semiclassical Condensation, PUSC, (Robnik 1988, Li and 
Robnik 1994, Robnik 1997), which claims the following: In the semiclassical limit 
^ — > the Wigner function of the n-th eigenstate condenses uniformly on the un- 
derlying classical invariant object (topologically transitive component) labeled by u, 
which can be an invariant N -torus, the entire energy surface (in case of ergodicity), 
or a chaotic component. The corresponding Wigner function is, in the most general 
case, given by (j^. This principle has a great predictive power, when accepted. 

The statistical properties of the wave functions in the coordinate space have been 
analyzed by Berry (1977b), where he has shown that in ergodic cases the proba- 
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bility amplitude distribution is Gaussian random functions, and its autocorrelation 
function is described by the Bessel functions. In case of integrability or quasi- 
integrability (a quantized invariant KAM torus) the wave functions in coordinate 
space are quite ordered, they typically have caustics (projection singularities when 
projecting the Wigner function (|1|) down the momentum space onto the configura- 
tion space). They are locally well described by the finite number of plane waves, 
because classically there is only a finite number of possible trajectory velocities (ob- 
tained by projection of the torus and of the quasiperiodic orbits on the torus onto 
the configuration space), whereas in case of ergodicity the number of plane waves 
is infinite (a circular ensemble of wave vectors), and they have uncorrelated phases, 
which implies Gaussian randomness. See e.g. (Robnik 1988, 1995, 1997). The limit- 
ing behaviour of the condensing Wigner function on a classical invariant object (no 
h enters in this equation, so in a sense we have classical Wigner functions!) implies 
also that the coarse grained probabihty density in the configuration space is just 
classical probability density, obtained by projecting the Wigner function onto the 
coordinate space, i.e. by integrating it over the momenta p. By coarse grained we 
mean smoothing over a few de Broglie wavelengths. And by the classical probabil- 
ity density we mean the value proportional to the time spent (asymptotically) in 
each cell of equal size (relative invariant measure) in the discretized phase space. 
Recently, we have brilliantly demonstrated this fact (Robnik et al 1998). 

^From the behaviour of the stationary wave functions we now turn back to the 
study of the statistical properties of the energy spectra, based on PUSC, giving 
some more details than in section 2. Again we restrict our discussion only to the 
time independent (autonomous) Hamiltonian systems with finite (bounded) classical 
motion and correspondingly a purely discrete quantal energy spectrum (no scatter- 
ing states). Further we asssume that there are infinitely many energy levels, so that 
the questions of statistical properties of energy spectra can be raised. Van Kam- 
pen (1985) has defined quantum chaos as " ...that property that causes a quantum 
system to behave statistically" . Now we have seen that this element is involved in 
the morphology of the eigenstates of classically chaotic, especially ergodic, systems. 
Therefore we must conclude that stationary quantum chaos exists, and corresponds 
exactly to the classical chaos. We shall now demonstrate that this is the case also 
when studying the energy spectra (and possibly also the statistical properties of 
other observables) . 

One of the most important cornerstones of the stationary quantum chaos is the 
so-called Bohigas-Giannoni-Schmit Conjecture (1984), BGS-Conjecture, introduced 
already in subsection 2.2, equation (p. It states that the classically ergodic Hamilto- 
nian systems (with discrete spectrum) exhibit universal spectral fluctuations, whose 
statistical properties are correctly captured by the conventional Random Matrix 
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Theories (RMT) (Mehta 1991), and are thus universal. See case (E) of section 2. 
If we ignore spin (which is not important in studying the classical limit) then the 
spectral fluctuations are described by either the fluctuations of the eigenvalues in 
the Gaussian Orthogonal Ensemble (GOE) of random matrices, if the system has 
an antiunitary symmetry, or by Gaussian Unitary Ensemble (GUE), if there is no 
antiunitary symmetry, such as e.g. time reversal symmetry, involved in the system. 
Of course, this statement applies to the statistical analysis of the spectrum, after 
the unfolding procedure, in which the actual physical energy is replaced by the av- 
erage number of states up to the given energy - i.e. the integrated (cumulative) 
level distribution. After unfolding the mean energy level spacing is by construction 
equal to unity. The average density of states typically is very well described by 
the familiar Thomas- Fermi rule of filling the classical phase space volume with the 
quantum cells of size (27r^)^, for which we have seen the reason in the above analy- 
sis of the semiclassical behaviour of the condensed Wigner functions of eigenstates. 
Sometimes the corrections to this asymptotically exact rule can be obtained, e.g. in 
plane and A^-dimensional billiards, constituting the famous Weyl rule with perime- 
ter, curvature, corner corrections etc (See e.g. Berry and Howls 1994 and references 
therein). In this unfolding procedure the information on the (nonuniversal) average 
density of states is eliminated from the spectrum, giving way to the possibility of 
universal fluctuations of energy levels around its nonuniversal mean distribution. In 
the classically ergodic systems this is exactly what we find, confirmed and supported 
by many numerical and actual experiments, and theoretically flrst corroborated by 
the result of Berry (1985) on the delta (A) statistics, and recently claimed to be 
proven by Andreev et al (1996), however, under much stronger conditions than er- 
godicity, namely assuming the exponential decay of correlations. For some recent 
review see (Robnik 1994,1995,1997). An important recent work in this connection 
is by Keating et al (1996,1997). 

What we said in the above paragraph applies to the scaling systems, where the 
energy limit (of quantum number n — > oo) is somehow equivalent to the semiclas- 
sical hmit of — > 0. That is, there is a scaling variable involving the energy and 
h such that the classical dynamics is constant while energy is changing. One such 
example are the billiard systems, among which the plane billiards are most widely 
used models. Examples of rigorously ergodic systems are the Sinai billiard, the sta- 
dium billiard of Bunimovich and the cardioid biUiard. The latter is the limiting case 
of the family of billiards with analytic boundaries defined as the quadratic confor- 
mal map of the unit disc, introduced by Robnik (1983,1984) and further studied by 
many workers. See Robnik et al (1997). In billiards the topology and the geometry 
of the phase portrait are exactly identical at all energies except for the scaling of 
the momentum as a square root of the energy. Therefore the limit of infinite energy 
is equivalent to the semiclassical limit oi k — > 0. The constancy of the classical 



16 



dynamics across the spectral stretches that we study is important to draw clear and 
safe conclusions about the relationship between the classical and quantum chaos. If 
a system is not a scaling system, then there is no way out other than taking just a 
small energy interval and letting — > so that in this limit the interval is containing 
an arbitrarily large number of energy levels, a necessary condition to introduce and 
to define statistical distributions. This general case is the one that we assumed in 
section 2. 

There are two most important statistical measures used to characterize the en- 
ergy spectra, and both of them are easily related to the E(/c, L) statistics that we 
introduced in section 2. One is the level spacing distribution usually denoted P{S), 
where ^S" is the length of the spacing and P{S)dS is the probability that S lies within 
the infinitesimal interval [S, S + dS) . It is normalized to unit probability (by defi- 
nition of probability density) and to unit first moment (due to the construction by 
the unfolding procedure). It can be shown that P{S) is the second derivative of 
the so-called gap probability E(0, L) of having no levels inside the interval of length 
S — L, so P{S) = d'^E{0, S)/dS'^. P{S) measures the short range correlations be- 
tween energy levels. Here the important point is the behaviour of P{S) at small S: 
For GOE one has the linear behaviour P{S) ~ const x S whilst for GUE we have 
P{S) const X S^. Correspondingly we talk about the linear and quadratic level 
repulsion: Because P{S) — > as 5" — > the level crossings, i.e. the degeneracies, 
are not likely, and in GUE this level repulsion is stronger than in GOE, giving rise - 
paradoxically - to a more regular spectrum. For GOE/GUE the corresponding P{S) 
for the infinite dimensional case cannot be obtained in a closed form, and for the de- 
tails the reader is referred e.g. to (Bohigas 1991). However, it is quite surprising and 
fortunate that the 2-dim GOE/GUE models yield closed analytic formulae which 
give an excellent approximation to the infinite dimensional case (which normally 
is referred to when speaking about RMT results). For GOE we have the so-called 
Wigner distribution 

PGOE(^) = ^exp(-^), (26) 

and for the GUE case 

on o2 A C2 

PGUE(-S) = ^exp( ). (27) 

71 TV 

Both cases are easily derived by assuming the 2-dim Gaussian real symmetric matri- 
ces (for GOE) or Hermitian symmetric matrices (for GUE). Interestingly, they can 
be derived from the so-called Wigner surmise (see e.g. Bohigas and Giannoni 1984, 
Brody et al 1981), which is an approximate argument outside the scope of RMT; it 
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is some kind of a statistical argument. 

Long range correlations are measured by the second most important statistical mea- 
sure, the delta statistics A(L), introduced by Dyson and Mehta. It is an inverse 
measure of spectral rigidity /regularity, as is immediately obvious from the definition: 

A(L) = {miuA BT / [-^{x) -Ax- Bfdx) (28) 

' L J-L/2 

where H{x) is the unfolded cumulative spectral staircase function {H{x) — linear 
average x plus oscillatory part M{x)), the minimum is taken with respect to the 
parameters A and and the average denoted by (...) is taken over a suitable en- 
ergy interval over x. Thus from this very definition A(L) is the local average least 
square deviation of the spectral staircase M{x) from the best fitting straight line 
over an energy range x of L mean level spacings. The more regular the (unfolded) 
spectrum, the easier is to find a linear fit over L levels, and consequently the smaller 
is A(L). The fact that we try to find the best linear fit to the spectral staircase 
implies that at small L the delta statistics A(L) always behaves as L/15 and thus 
carries no information about the system at all. 

Sometimes it is useful also to know the number variance, denoted by Y?{L), the 
dispersion of the number of levels n{L) in an interval of length L, where {n{L)) — L, 
and 

E2(L) = ((n(L) - L)2) = L - 2 [\l - r)Y2{r)dr (29) 

Jo 

where Y2{r) is the pair cluster function (Bohigas 1991, Haake 1991, Mehta 1991). 
There exists also the connection between and A, 

A(L) = ^ ["^{L^ - 2LV + r^)J:\r)dr (30) 
Jo 

although, strictly speaking, this has been proven so far only within the context of 
RMT (Aurich et al 1997, Mehta 1991). Finally, as announced in section 2, we shall 
consider not the set of all the cluster functions Yn{xi, 2:2, ... , Xn), where n = 2, 3, . . ., 
but rather the E{k, L) statistics, for all A; = 0, 1, 2, . . ., following the suggestion of 
Steiner and coworkers (Aurich et al 1997), because they are very easy to calcTilate 
numerically and yet contain the complete information about the spectral statistics. 
Since by definition E(A;, L) is the probability that inside an interval of length L we 
find exactly k levels, there are simple relationships to other statistical quantities. 
For example, as already mentioned, the level spacing distribution P{S) is 
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P{S) 



E(fc = 0,L = ^) 



(31) 



and 



T?{L) = Y.{k-Lfnk.L). 

k=0 



(32) 



and therefore, through (|30|) , we have the relation expressing A(L) in terms of the 
E(/c, L) statistics. 

/^From the definition of the Poissonian E(fc, L) statistics in (0) it is easily derived 
(see m) 



and after (^) 

and then using (RDI) 



PpoissoniS) = exp(-S') 



(33) 
(34) 



^Poissoni^L) . 
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(35) 



Poissonian statistics means also by definition that there are no correlations, i.e. the 
pair correlation function factorizes, so that we have for the pair cluster function (c.f. 
Mehta 1991, Bohigas 1991) 



Poisson 



(x) = 0. 



(36) 



Thus, using this fact in equation (|29|) and then (|30|) we again recover Poissonian 
values (|3^) and (^). 

It is at large L that the different universality classes of behaviour emerge. Interest- 
ing are the asymptotic results for large L. For a completely regular (equidistant) 
spectrum (like e.g. one-dimensional harmonic oscillator) one obtains that for large 
L ^ 1 , A(L) is just constant and equal to 1/12. On the other hand the RMT gives 
for GOE 



Agoe(^) ~ ^logL, 



(37) 
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and for GUE 



AGVE{L)^—\ogL, (38) 

This result of RMT has been derived on dynamical grounds for the energy spectra of 
individual classically ergodic systems, applying the Gutzwiller's periodic orbit the- 
ory, in a remarkable and important paper by Berry (1985), giving some theoretical 
support to BGS-Conjecture. As it is believed that the level repulsion is a purely 
quantal effect and cannot be derived semiclassically (Robnik 1986, 1989), one is not 
surprised that for eleven years there was hardly any theoretical progress in estab- 
hshing the BGS-Conjecture. Indeed, Berry concluded (1991) that P{S) cannot be 
derived from applying periodic orbit theory, because it depends sensitively on orbits 
of all lengths (periods). However, recently Andreev et al (1996) claim to have derived 
BGS-Conjecture by a different thinking, namely by studying the spectrum of the 
Frobenius-Peron operator in the semiclassical limit, using the techniques from the 
supersymmetry field theories, especially the nonlinear sigma model (Weidenmiiller 
et al 1985). In considering the semiclassical limit they assume stronger properties 
than ergodicity, namely the exponential decay of (classical) correlations. 

We may conclude that the support for the BGS-Conjecture is so strong that it 
can be regarded as well established, although not rigorously proven as yet. By this 
I mean especially the unusually strong and massive numerical support and evidence 
accumulated during the past fourteen years. Thus given the correctness of BGS- 
Conjecture we speak about the universality classes of spectral fluctuations, namely 
the GOE and GUE class, of subsection 2.2. As is seen in the above formulae for 
P{S) and for A(L), and in the general equation (^, the universality is indeed estab- 
lished: There is no parameter in the statistical properties of spectral fluctuations, 
and the statistical measures are identical for all ergodic systems, irrespective of their 
dynamical and geometrical details. Turning this aspect around, we conclude that 
the spectral fluctuations in a classically ergodic system do not have any further in- 
formation content. When talking about BGS-Conjecture as applied to individual 
classically ergodic systems we must emphasize that if there are any exact unitary 
symmetries involved in the system, then they must be first eliminated before ap- 
plying BGS-statement. This process we call desymmetrization. Then, even after 
desymmetrization, we still have to decide whether GOE or GUE statistics apply: 
The general classification criterion is: If the system has an antiunitary symmetry 
(like e.g. the time reversal symmetry) then GOE applies, and GUE otherwise (i.e. 
if there is no antiunitary symmetry) (Robnik and Berry 1986, Robnik 1986). 

It is a surprise that RMT apply so well to individual dynamical systems. We claimed 
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that BGS-Conjecture holds true in the strict semiclassical hmit ^ — ^ 0. However, 
how small must be h to see this happening? Thus we are now addressing the ques- 
tion of the limitations to universality due to the not-sufficiently-small value of the 
(effective) Planck constant. 



The following criterion is important: As soon as the eigenstates are fully extended 
chaotic (ergodic) in the sense of (Q) the BGS-Conjecture applies. This is always 
happening asymptotically, as ^ — >• 0, but a rough criterion is that the classical diffu- 
sion time (the typical time for the classical dynamics to conquer the entire available 
phase space - energy surface) is shorter than the break time introduced in section 
2.4, equation (12): If this inequality applies (strong enough, i.e. by a factor 10 or 
so) then the eigenstates will be fully extended. If the inequality is reversed, then the 
eigenstates are chaotic (they lie in a classically chaotic region) but localized, which 
means occupying only a small piece (a proper subset) of the dynamically available 
phase space. Obviously the break time {27rh) / (AE) goes to infinity a.sh—>-0, whilst 
the classical transport time is independent of h, and thus the desired inequality is 
asymptotically always satisfied, and therefore in the limit all semiclassical states are 
fully extended states and we recover the universality of section 2. 



In such a dynamically localized classically ergodic regime another interesting phe- 
nomenon occurs, namely the fractional power law level repulsion, by which we mean 
that 



P{S) ^ const X S^, (39) 

where the exponent (3 can be 0, 1, 2 or anything in between. Thus the localization 
phenomena of the chaotic eigenf unctions soften the strength of the level repulsion 
{(3 going from 1 to 0). This phenomenon seems quite obvious, since the tails of the 
localized wave functions overlap even less and thus interact less strongly resulting 
in reducing the value of (3. This trend towards the Poisson statistics is well known 
in the context of localization phenomena in disordered solid state systems. At 
present we do not have a theory on how f3 should be related to the localization 
lengths/areas/volumes, and how to predict them. But we have numerical examples 
demonstrating these features (see e.g. Prosen and Robnik 1994a,b). The theory 
must satisfy the known limit of — 1 when h ^ 0. One distribution function 
which at present does not have any deep physical justification as yet, but is just a 
nice mathematical model which captures the global level spacing distribution with 
the local property (|39| ) is the well known Brody distribution (Brody 1973, Brody et 
al 1981) 



21 



PBroAj{S;f3)=aS^exp{-bSP+'), a = b{(3 + l), b = {r{^^)y+' (40) 

where a and b are obviously determined by the normahzations of the total probability 
and of the first moment to unity. For jS = we have Poisson distribution (expo- 
nential, i.e. P{S) = exp(— S*)), and for /3 = 1 we have Wigner (i.e. 2-dimensional 
GOE, given in equation (|26|)). The role of dynamical localization phenomena has 
been first realized by Chirikov et al (1981) in the time-dependent systems like kicked 
rotator and Rydberg atoms in microwaves, but has been also suggested in the time- 
independent Hamiltonian systems (Chirikov 1993). Feingold (1996) has recently 
found deeper relationship between the two phenomena. 

After having explained the two universality classes of spectral fiuctuations in the 
classically ergodic systems, we now have to add and rediscuss the third universal 
class comprising of classically integrable systems of two or more degrees of freedom 
(Robnik and Veble 1998), see subsection 2. (Systems with only one degree of free- 
dom are exceptional and special in the sense that as /i — > the local spectrum is 
just the perfectly regular equidistant spectrum.) Indeed, if we have two or more 
quantum numbers the entire spectrum can be thought of as being composed of an 
infinite number of statistically uncorrelated number sequences, which of course must 
result in the Poisson statistics E(A;, L) in (0), and specifically (|33|), (|3^ ) and (pS]). 

There are semiclassical arguments resting upon the torus quantization by Berry 
and Tabor (1977) showing that the statistics should be Poissonian. The torus quan- 
tization (EBK quantization) is embodied in equation (|lj), describing the associated 
Wigner functions. However, these semiclassical arguments are only approximation 
and it is far from obvious that they should correctly describe e.g. the fine structure 
of energy spectra and thus the level repulsion and their absence, as has been recently 
pointed out by Prosen and Robnik (1993c, d). In fact, as explained above, it is be- 
lieved that semiclassics cannot explain the level repulsion (short range correlations), 
since it is a purely quantum effect. We know that there are exceptions from the 
Poissonian behaviour which have been rigorously proven to exist and involve some 
highly nontrivial and sophisticated mathematics (Bleher et al 1993). Nevertheless, 
there is quite massive numerical and experimental support to the statement that the 
spectral fiuctuations of classically integrable systems are quite accurately described 
by the Poisson statistics (^ (Robnik and Veble 1998). There might be cases where 
the statement is rigorous, whereas in general we think that the measure of excep- 
tions is small and maybe vanishing in some sense. In every case this delicate problem 
persists to be very important and difficult, but it should also be accepted that typ- 
ically Poisson model is an excellent approximation. The most important feature 
is the absence of short range correlations implying the absence of level repulsion. 
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which means that degeneracies are allowed and this is mathematically exhibited in 
P{S) —>■ const 7^ 0, in fact according to (RSp we have P{S) ^ 1 as 5* ^ 0. 



Another remark should be made about a limitation to universality, mentioned in 
subsection 2.4, the behaviour of A(L) in individual dynamical systems at large L, 
where the limitations to universality of section 2.4 set in. As discovered by Casati 
et al (1985) and later explained by Berry (1985) there is the phenomenon of satu- 
ration, by which we mean that A(L) becomes effectively constant and equal Aqo if 
L > Lmax, in O'l^y system (ergodic, integrable, partially chaotic - KAM type). This 
leveling off of the delta statistics is nonuniversal, but the Lmax can be estimated in 
the context of Berry's theory (1985) as in equation (|lTl). Therefore for any N > 2 
the onset of saturation Lmax goes to infinity in the semiclassical limit h 0, giving 



way to the full universality of the three universality classes (|37| , pSj |35|) . The details 



of the saturation value A^o at a fixed and nonzero h can be found in (Berry 1985). 

Finally, we should say something more about the P{S) and A(L) in the classi- 
cally mixed systems, thus giving more details of the Berry- Robnik (1984) picture. 
Both statistics are of course implied by the most general statistics (^, through the 
formulae (|31|) for P{S), and through (|3^) and ( pOD for sigma and delta statistics. 



Using the general equation (H) and approximation (p6|) to first calculate Egoe, 
and then to find P{S) through (|3l|) , for m level sequences, we obtain an explicit 
analytic formula for P{S), namely 



j2 m rz 

Pm{S) = Js-2^^M-PiS) n erfc(^p,5)] (41) 

where erfc(x) = J!^ dt expl—t"^) is the complementary error function. (For GUE 
one has to use (|27|) instead of (EB).) One can show 



PmiS = 0) = l-J2pl (42) 
i=2 

so that now as a consequence of the statistically indpendent superposition of a num- 
ber of level (sub) sequences we get a trend towards Poissonian statistics, since the 
degeneracies become possible due to the lack of level repulsion among the regular 
levels on the one hand and among the levels belonging to different sequences. Thus, 
e.g. even if there is no regular component but two equally strong chaotic compo- 
nents, so that pi = 0, but p2 = Ps = 1/2, we get PmiS = 0) = 1 - 1/4 = 3/4 ^ 0. 
Only if there is only one chaotic (ergodic) component we find Pm{S = 0) = 0, de- 
scribing the GOE-like level repulsion. 
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Most important in practical applications is the 2-component Berry-Robnik formula 
(m = 2), namely 



P2{S,pi) = plexp{-piS)eTic{^p2S) + (2pip2 + ^TTplS) exp{-piS - ^TrplS^), 

(43) 

and we see the special case of equation (|4^ ) 

P2{S = 0,p^) = l~pl = p^i2-p^), (44) 

which vanishes only iff pi = and p2 = 1 (ergodicity). This level spacing distribu- 
tion is very important especially in practical applications, because in mixed systems 
typically we have a very large dominant chaotic region, so that the next largest 
chaotic region is much smaller by orders of magnitude, say only one percent of the 
leading one and can be neglected. In such case the two component formula (^) is 
an excellent approximation. 

For the delta statistics (^) one can derive the additivity property implied by the 
statistical independence of the (sub) sequences. First one shows it for the sigma 
statistics (the number variance) (see e.g. Bohigas 1984,1991). Then one can show 
(Seligman and Verbaarschot 1985) 

m 

A(L) = 5:A,(p,L). (45) 
i=i 

It is now interesting to verify whether Berry-Robnik theory applies to actual sys- 
tems which we can analyze numerically. We (Prosen and Robnik 1993c, 1994a,b) 
have done such analysis for a certain one-parameter family of billiards, introduced 
by Robnik (1983, 1984), namely the 2-dim billiard shape defined as the complex 
quadratic conformal mapping w = w{z) = z + Xz"^ of the unit disc \z\ < 1 in the 
2;-plane onto the w-plane. The boundary curve is known in the theory of analytic 
curves as the Pascal's Snail. (Usually workers refer to this billiard system as Rob- 
nik billiard, because it was introduced and dynamically analyzed by the author.) 
The system is very important because it has analytic boundaries up to the limiting 
value of the shape parameter A = 1/2 where the singularity appears aX z = —1, and 
therefore for small A the KAM-Theory applies: At A = we have the integrable 
case of the circle billiard with conserved angular momentum. For small A < 1/4 we 
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have a convex billiard with analytic boundaries, studied extensively especially by 

Lazutkin (1981,1991), where much can be said about the classical and semiclassi- 
cal analysis, including a construction of an approximate integral of motion (Robnik 
and Berry 1985, Robnik 1986). This is essentially KAM scenario. For A > 1/4 
the boundary is nonconvex and the KAM theory does not apply because the bounce 
map becomes discontinuous. When A = 1/4 the first point of zero curvature appears 
at z = —1 and according to Mather (1982) this guarantees that all the Lazutkin 
caustics (generated by invariant tori for glancing orbits supporting the whispering 
gallery modes of quantal eigenstates in the semiclassical picture) are destroyed, giv- 
ing way (preparing the way) for ergodicity, which has been postulated by Robnik 
(1983). In fact, a careful analysis of certain periodic orbits (Hayh et al 1987) has 
shown that also for A > 1/4 they can be stable, surrounded by very tiny stability 
islands that can hardly be detected numerically, which was the reason why in the 
early work (Robnik 1983) they have not been seen. According to their estimates 
the system has stable islands up to A 0.2791, whilst Li and Robnik (1996) have 
numerical evidence that ergodicity is possible for A > 0.2775. Recently it has been 
rigorously proven by Markarian (1993) that for A = 1/2 (we have a cusp singularity 
at ^ = —1, because dw/dz = there) the system is ergodic, mixing and K. This is 
thus the first billiard system with analytic boundaries having the chance to be er- 
godic for 0.2775 < A < 1/2. (The Sinai biUiard and the stadium of Bunimovich are 
rigorously ergodic, mixing and K, but they do not have analytic boundaries.) The 
system has been recently studied by many workers (Berry and Robnik 1986, Robnik 
and Berry 1986, Frisk 1990, Bruus and Stone 1994, Stone and Bruus 1993a,b, Backer 
et al 1995, Bruus and Whelan 1996). 

The system is thus ideal to study the morphology of quantum eigenstates at various 
A, following a continuous transition from the domain of torus states and Poisson 
statistics, through the regime of the generic behaviour of mixed dynamics, to the 
extreme case of (rigorous) ergodicity and entirely chaotic states with GOE statistics. 
(Of course, as explained in the introduction, we must separate the exact symmetry 
classes of a given dynamical system before performing the statistical analysis of the 
energy spectra. This procedure is called the desymmetrization. In our case we have 
even and odd reflection symmetry classes.) The main results on this have been 
published in (Li and Robnik 1995a,b,c) 

For the details please see (Robnik 1997) and the references therein. Another bril- 
liant numerical conflrmation of the Berry-Robnik statistics was given in (Prosen 
and Robnik 1994a, b) for the quantized compactifled standard map, and by Prosen 
(1995,1996) for a semiseparable oscillator, and recently for a quartic billiard (Prosen 
1998). The main problems and issues cocerning the Berry-Robnik picture have been 
expounded recently in our comment (Robnik and Prosen 1997). The deviation from 
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Berry-Robnik regime towards the Brody-like behaviour with fractional power law 
level repulsion, as described in subsection 2.4, has been analyzed in detail in (Prosen 
and Robnik 1993c, 1994a,b). 



4 Statistical properties of classically chaotic mo- 
tion and the measure of chaotic components 

In this section we address the question of the statistics of classically chaotic mo- 
tion and the problem of how to determine the measure of the chaotic components 
P2,p3, ■ ■ which in this section we shall simply denote by p2, relevant for the prob- 
lem of stationary quantum chaos in mixed systems, dealt with in subsection 2.3. 

In a recent work (Robnik et al 1997) we have demonstrated some general scaling 
laws in the behaviour of stochastic diffusion in strongly chaotic systems (ergodic, 
mixing and K with large Lyapunov coefficient, i.e. large KS entropy), mainly in 
Hamiltonian systems, or in the strange attractors of dissipative systems. The so- 
called random model that we developed describes very well the diffusion on chaotic 
components, in the sense that the relative (invariant) measure p{j) as a function of 
the discrete time|^ j approaches unity exponentially as 



where A''^ is the number of cells of equal size (relative invariant measure) a = l/N^ 
into which the whole ergodic component is decomposed, provided A*";, is sufficiently 
large, say Nc > 100 or so. In the above equation we have defined p(j) = P2(j) / P2(oo) . 
Thus the average measure of occupied domain on the grid of cells is (ka) = p{j). 
This random model rests upon the assumption that there are absolutely no correla- 
tions, not even between two consecutive steps, so that at each step (of filling the A^^^ 
cells) we have the same a priori probability a = 1/Nc oi visiting any of the cells, ir- 
respective of whether they are already occupied or not. Such absence of correlations 
can be implied and expected by the large Lyapunov exponents, which in turn imply 
strong stretching and folding (of a phase space element) even after one iteration, 
meaning that such a phase element will be evenly distributed (in the coarse grained 
sense) over the entire phase space (or surface of section). The universal scaling 
property is refiected in the fact that p(j) is a function of the ratio (j/Nc) only, and 
does not depend on j and Nc separately. 

^We work either with mappings or with Poincare mappings on the surface of section. In each 
case j is the number of the iterations of the map. 




(46) 



26 



Such an assumption of absence of all correlations appears to be strong at first sight, 
and therefore it is quite surprising that the model describes a whole lot of deter- 
ministic dynamical systems for which we can expect large Lyapunov exponents, 
namely 2D bilhard (Robnik 1983, A = 0.375), 3D bilhard (Prosen 1997a,b, a = -i, 
b = — y), ergodic logistic map (tent map), hydrogen atom in a strong magnetic 
field (e = —0.05) (Robnik 1981, 1982, Hasegawa et al 1989), and standard map at 
[k = 400), in which the agreement is almost perfect, except for the last two systems 
where we see some long-time deviations on very small scales. However, in the stan- 
dard map at k = 3, and in Henon-Heiles (1964) system at E = ^ the deviations are 
noticeable though not very big (about only 1%). 

It is also quite astonishing that the random model applies very well even to ergodic- 
only systems, with strictly zero Lyapunov exponents, namely in case of the rectangle 
billiards (Artuso et al 1997), where the deviations from the exponential law (^) on 
the largest scale is within a few percent only. It is a well known result (Sinai 1976) 
that polygonal billiards have exactly zero Lypunov exponents, easy to understand 
since all periodic orbits are marginally stable (parabolic), and since they are every- 
where dense, we conclude that the Lyapunov exponents must be zero everywhere. 

As a small but interesting comment we should mention our results on testing the 
random number generators from (Press et al 1986), where two of them (ranO and 
ran 3) are found to be in perfect agreement with the random model, whilst the other 
two (rani and ran2) are exhibiting big deviations. Thus, indeed, some deterministic 
dynamical systems like hydrogen atom in strong magnetic field etc. can be better 
number generators than some built-in (black-box) computer algorithms. It should 
be acknowledged, however, that there are other random number generators which 
pass all tests of randomness, including ours, e.g. in (Finocchiaro et al 1993). 

The random model developed in (Robnik et al 1997) is a statistical model which 
predicts not only the average relative measure of occupied cells p(j) = {ka), the 
average taken over k, in (P6|), but also the standard deviation cr(j), which under the 
same assumption of sufficiently large Nc is equal to, to the leading order. 



^(J) = \/((H^)-((H)^ = \/^^, (47) 
and gives us an estimate of the size of expected statistical fluctuations in p(j). 

The random model (Robnik et al 1997) has been subsequently generalized in an 
important direction (Prosen and Robnik 1998), namely to describe the diffusion on 
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chaotic components in systems of mixed dynamics, with divided phase space, having 
regular regions (invariant tori) coexisting in the phase space with chaotic regions, a 
typical KAM scenario (Kolmogorov 1954, Arnold 1963, Moser 1962, Benettin et al 
1984, Gutzwiller 1990). Such systems in two degrees of freedom can have the fractal 
boundary between the regular and irregular component and thus the convergence 
to the theoretically expected results can be very slow, mimicking a departure from 
the random model, although ultimately it conforms to this model. In three or more 
degrees of freedom there is no boundary between the regular and chaotic regions, 
because we have the Arnold web (Chirikov 1979), which is everywhere dense in the 
phase space, and thus a naive box-counting would imply always that the relative 
invariant measure of the chaotic component is equal to the measure of the entire 
phase space, so P2(j) = 1, which is wrong, because the KAM theorem gives rig- 
orously that the relative measure of the regular component pi is strictly positive, 
pi > 1, moreover it is close to unity with the perturbation parameter. We assume 
that the invariant measure of the chaotic component is positive, although strictly 
speaking this is a major open theoretical problem in the mathematics of nonlinear 
systems, the so-called coexistence problem (Strelcyn 1991). Therefore in such case 
one must introduce the possibility of different a priori probabilities, which now are 
no longer just the same and equal to a = 1/iVc, but have a certain distribution de- 
scribed by the so-called greyness distribution w{g), where 5^ is a continuous variable 
on the interval [0, 1]: g = means no visits (white cells), g = 1 are the most fre- 
quently occupied cells (black cells), and those cells with < g < 1 have intermediate 
number of visits (grey cells). With this model we have shown how by measuring 
(numerically calculating) w{g) we can determine the relative invariant measure p of 
the chaotic component. The result is 

p = / gw{g)dg (48) 
Jo 

and the time dependent relative measure of occupied domain is equal to 



P(j) = 1 - dgw{g) exp(-^) (49) 
JO pAc 

and the standard deviation is still given precisely by the equation (^). The greyness 
distribution can be calculated numerically quite easily by noticing that the greyness 
g is proportional to the average occupancy number n{g), namely n{g) = ^f/p, so by 
measuring n{g) in the limit j — * 00 and after normalizing the g of the peak of n{g) 
to unity, we get the g^s, and then by binning them into bins of suitably small size 
A^f we get the histogram for w{g). 
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In case of ergodicity (only one chaotic component) we have w{g) = 6{g — 1), the 
Dirac delta function at g = 1, and then from equations and (^) follows the 
random model, with exponential behavior (|46|). 



In a later work we dealt with only ergodic systems, but such having several compo- 
nents, each of them also being ergodic, but weakly coupled, by which we mean that 
the transition probability for going from one to another component is very small 
and the typical transition time j* very long. Obviously, at small times we shall find 
the random model (M) with being equal to the number of cells of the starting 
component, Nc = Ni, whilst for very large times j, bigger than the typical transition 
time j*, so j ^ j*, we shall find again the random model (|46|) , but now with Nc 
being equal to the number of all cells in the system, Nc = Ns- In between, when 
j ^ j*, we have the crossover regime which we analyse in the present work. So, the 
finite time structure of ergodic systems controlled by their transport times can be 
also captured analytically. We call the corresponding model the multi-component 
random model of diffusion. For more details see (Robnik et al 1998). 



As the concluding remark of this section it should be explained that the relative 
measures pi, P2, ■ ■ ■ , Pm we need in the formula for the statistics of the classically 
mixed systems in subsection 2.3, equation (^), are the relative Liouville measures on 
the energy surface (E), the surface being defined by H{q, p) = E = const., whereas 
in this section we explained how to calculate the relative (invariant simplectic) mea- 
sure /i of a given invariant chaotic set u on the Surface of Section (SOS). They 
are not the same, for a given chaotic set, but the procedures to calculate them are 
relatively simply related to each other, as has been explained by Meyer (1985). The 
answer is obtained from the general relation 



{A)e = I VpA(q, p)5,{E - ff(q, p)) = {tA)sos = f dXrAiq, p), (50) 
Je Jsos 

where dX is the invariant simplectic measure on SOS, saying, that the microcanon- 
ical average of any classical function A(q, p) (observable) over the energy surface E 
is equal to the average over the SOS of A weighted by the average time of recurrence 
r of a trajectory to SOS on each invariant (ergodic) component, r is constant for 
a given trajectory and thus also inside a given invariant (ergodic) component, but 
changes from one component to another, also on irrational invariant tori. Thus if, 
specifically, A^{c[, p) is the characteristic function on a chaotic invariant component, 
so equal to Xw(q, p), then we have for the measure p(u;) of equation (H) 
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Jj,d^qd^p6,{E-H{q,p)) 
Isos dXr Jsos d^r 



(51) 



The third equahty in the above equation follows by recalling that r is constant on 
a given invariant component u, equal to t^. Of course, we can and should calculate 
(numerically) p{uj) by discretizing the SOS into a network of Nc cells of equal size 
a = l/Nc- As is evident from equation (0), the invariant measure dX in a given 
cell can be chosen simply equal to 



This solves our problem of determining the classical measures p{uj) in the context of 
their role in the semiclassical limiting behaviour of the eigenstates and of the Hilbert 
space of a general, classically mixed, and therefore generic system, dealt with in the 
context of quantum chaos mainly in subsection 2.3. 



5 Discussion and conclusions 

The main purpose of this paper is first to provide a compact review of the main 
topics in the stationary quantum chaos in the general quantal Hamiltonian systems 
with discrete energy spectra, in correspondence with classical chaos, in the strict 
semiclassical limit where the effective Planck constant h is sufficiently small. This 
problem has been expounded in the Introduction, section 1. 

In section 2, we have explained the universality classes of spectral fluctuations and of 
their wave functions and of Wigner functions. These comprise of the classically in- 
tegrable systems exhibiting Poisson statistics (subsection 2.1), and of the classically 
ergodic systems, obeying the statistics of the of the eigenvalues of the ensembles 
of random matrices from the classical Random Matrix Theories, namely GOE and 
GUE, depending on whether the system has or not an antiunitary symmetry, like 
e.g. the time reversal symmetry (subsection 2.2). Then we went on explaining the 
general case of classically mixed dynamics, in the transition region between classical 
integrability and ergodicity, very often well described by the scenario of the KAM 
Theory (subsection 2.3). We have shown how the Berry- Robnik (1984) picture can 
be applied, and have generalized it to arbitrary statistics, speciflcally E(A;, L) statis- 
tics. The main message is that the Hilbert space of eigenstates of a mixed system is 



dX = gw{g)dg 



(52) 
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decomposed into the set of regular states (associated with a EBK/Maslov quantized 
invariant torus) and a set if irregular (chaotic) eigenstates, where their properties 
are captured by the most general semiclassical Wigner function in equation (|^). The 
regular and irregular sequences are split exactly in proportion to the relative invari- 
ant Liouville measures of their supporting classical invariant sets. Then we have 
described some limitations to the universality classes and the semiclassical asymp- 
totical behaviour (subsection 2.4). The first limitation stems from the existence of 
the Berry's outer energy scale (after unfolding) L^ax = {^nh) / (AE) (see equation 
(|TI|)), so that at L > L^ax we have no longer universality but saturation (e.g. of 
the sigma and of delta statistics). Nevertheless, as ^ ^ 0, Lmax goes to infinity. 
The second limitation comes from the localization phenomena, controlled by the 
relation of the two time scales, the break time and the classical diffusion time. If 
break time for a given stationary eiegnstate is shorter than diffusion (ergodic) time, 
then we have a strong localization. If the system is ergodic (but slowly diffusing), 
then we find a departure from the universal RMT statistics, and in extreme case 
can be close to Poisson statistics. In the opposite extreme, we find the extended 
states and correct behaviour generally described by the mixed case (M) of subsection 
2.3. In subsection 2.5 we give some fundamental results on the statistical properties 
of general matrix elements, generalizing the important Feingold-Peres (1986) theory. 

In section 3 we discuss the fundamental propeties of the Wigner functions of the 
eigenstates, and introduce the Principle of Uniform Semiclassical Condensation. 
Then we say more about the wave functions and the spectral statistics, especially 
about the level spacing distributions, sigma and delta statistics. There we have also 
recalled Van Kampen's definition of quantum chaos and concluded that the quantum 
chaos, as a phenomenon " ...causing the quantum systems to behave statistically..." , 
does exist in the problem of the stationary Schrodinger equation and its solutions. 

In section 4 we deal with the problem of how to determine, also numerically, the 
relative invariant measure of classically chaotic components, which enter in the semi- 
classical formulae for statistics of mixed systems (^. We have described the main 
ideas, the approach and the results of our recent works on this subject, giving the 
final decription of how to proceed. The main goal of the theory of stationary quan- 
tum chaos is to explain and to describe the mixed systems. Then, the quantal 
parameters pj determined by analyzing the spectral statistics must be equal to the 
purely classical relative invariant (Liouville) measures (subsection 2.3). 

We believe that the material presented here is also a stimulation for further the- 
oretical and numerical work, and also shows the need for more rigorous results on 
the side of the mathematical physics. For example, we still need the proof of the 
BGS-Conjecture, especially in its full generality presented in equation (^. We need 
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to prove the Principle of Uniform Semiclassical Condensation, which is more general 

than the BGS-Conjccture. Finally, as one of the most important issues, we need 
to understand more deeply and quantitatively the localization phenomena in the 
stationary eigenstates of general systems. Some important new steps in this direc- 
tion have been recently undertaken by Casati and Prosen (1998a,b), and by Krylov 
and Robnik (1998). The most general aspects of quantum chaos related to other 
branches of theoretical and experimental physics have been recently reviewed and 
discussed in detail by Weidenmiiller and coworkers (Guhr et al 1998). 
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